in 


ABERDEEN  PROVING  GROUND,  MARYLAND 


ASTIA  AVAILABILITY  NOTICE 

Qualified  requestors  my  obtain  copies  of  this  report  from  ASTIA. 


The  findings  In  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Anqr  position. 


BALLISTIC  RESEARCH  LABORATORIES 


REPORT  NO.  1191*- 
JANUARY  1963 


A  SOLUTION  OF  THE  GODDARD  IROBIBM 


Boris  Garfinkel 


Computing  Laboratory 


RDT  &  E  Project  No.  1M010501A003 


ABERDEEN  PROVING  GROUND,  MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


REPORT  NO.  1194 


BGarfinkel/bj 

Aberdeen  Proving  Ground,  Ml. 
January  1963 


A  SOLUTION  OP  THE  GODDARD  IROBLEM 

ABSTRACT 

The  problem  of  optimizing  the  thrust  of  a  vertically  ascending  rocket 
is  solved  here  under  the  assumption  of  isothermal  atmosphere  in  two  important 
cases:  1)  the  jet  Mach  number  is  sufficiently  large;  2)  the  drag  is  a 
convex  function  of  the  velocity. 

The  first  case  embraces  all  physical  drags  and  is  valid  for  the  Earth; 
the  second  extends  to  all  atmospheres,  but  is  restricted  to  drags  that  are 
fairly  common. 

With  impulsive  boosts  in  velocity  admitted,  the  solution  is  shown  to 
contain  a  finite  number  of  such  boosts  in  the  sonic  region  of  the  rocket 
velocity,  and  to  contain  no  coasting  arcs  except  in  the  terminal  stage. 

An  absolute  minimum  is  proved  with  the  aid  of  a  Sufficient  Condition 
applicable  to  problems  of  optimum  contrpl. 


3 


TABLE  OP  CONTENTS 


Page 

1.  INTRODUCTION  .  7 

2.  FORMULATION  OF  THE  PROBLEM .  7 

3.  THE  AUXILIARY  PROBLEM .  9 

1*.  SOME  PROPERTIES  OF  THE  DRAG .  10 

5.  THE  EULER  EQUATIONS .  13 

6.  THE  TRANSVERSALITY  CONDITION .  15 

7.  THE  CORNER  CONDITION .  16 

8.  THE  HILBERT  CONDITION .  18 

9.  CONDITIONS  OF  LEGENDRE  AND  WEIERSTRASS .  19 

10.  THE  JACOBI  CONDITION .  20 

11.  THE  SUFFICIENCY  CONDITION .  22 

12.  THE  STEADY  STATES  OF  MOTION .  22 

13.  THE  BASIC  THEOREMS . 24 

14.  SUMMARY .  26 

APPENDIX .  27 

REFERENCES . .  .  28 


5 


1.  INTRODUCTION 


The  problem  of  maximizing  the  summit  altitude  of  a  vertically  ascending 
rocket,  of  which  the  Goddard  Problem  (1919)  is  a  variant,  has  received 
considerable  attention  in  the  literature.  One  of  the  earliest  attacks  on  the 
problem  should  be  credited  to  Lewy  (1944),  instigated  by  Dr.  R,  H.  Kent,  then 
the  Associate  Director  of  the  Ballistic  Research  Laboratories.  Despite  the 
notable  advance  achieved  by  Tsien  and  Evans  (1951) >  numerous  gaps  in  the  theory 
still  remain  to  be  filled.  As  has  been  pointed  out  by  Leitman,  Ross,  et  al., 
the  problem  continues  to  be  beset  by  the  difficulty  arising  from  the  require¬ 
ment  that  the  mass  be  monotone.  Solutions  that  meet  this  requirement  have 
been  obtained  only  in  a  few  very  special  cases,  typified  by  the  work  of 
Miele  (1958),  who  treated  flight  in  vacuum  and  the  power  law  of  drag. 

In  the  present  paper,  which  is  an  outgrowth  of  the  author's  unpublished 
work  of  1949,  reported  at  a  Ballistic  Research  Laboratories  Colloquium,  the 
class  of  soluble  cases  is  considerably  broadened.  With  the  assumption  of 
isothermal  atmosphere  and  the  admissibility  of  infinite  thrust,  a  solution  is 
obtained  in  the  following  two  cases: 

1)  The  jet  Mach  number  is  sufficiently  large. 

2)  The  drag  is  a  convex  function  of  the  velocity. 

The  first  case  is  valid  for  the  Earth;  the  second  is  restricted  to  a  class  of 
drags  that  are  fairly  common.  The  remaining  case,  where  neither  (l)  nor  (2) 
holds,  is  being  left  as  a  subject  for  future  investigation. 

A  recapitulation  of  the  relevant  existing  theory,  designed  to  provide  the 
necessary  background  for  the  current  development,  is  incorporated  in  sections 
2  and  5. 


2,  FORMUIATION  OF  THE  PROBLEM 

The  equation  of  motion  of  the  rocket,  subject  to  forces  of  gravity,  drag, 
and  thrust,  is 

me  +  mV  +  |  CD(V,X)p(X)V2S  +  mg(X)  =  0,  (l) 
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where  m  is  the  mass,  V  the  velocity  of  the  rocket,  the  drag-coefficient, 

X  the  altitude,  p  the  density  of  the  air,  S  the  cross-section,  g  the  acceleration 
of  gravity,  c  the  jet  velocity,  and  the  superscript  dot  indicates  the 
differentiation  with  respect  to  the  time. 

We  shall  introduce  the  simplifying  assumptions: 


1)  CD  is  a  function  of  V  only, 

2)  g  =  const. , 

3)  P  =  PQ  exp  (-x/i),  t  =  const., 

define  the  dimensionless  parameters  a,  £3  by 

2  ,  2 
a  =  gi/c  ,  P  h  2  m08/c  PQS, 

0  <  a  <  oo  ,  0  <  (3  <  oo  , 

where  m  is  the  initial  mass,  and  dimensionless  variables  x,  v,  a>, 
o 

X  =  gX/c~  ,  V  =  v/c  ,  CD  =  log  mo/m  , 
y  ~  cd  —  v  —  x/a  ,  f  =  CDv  eV/P  • 

Then  (l)  becomes 


(2) 


(5) 

y,  and  f  by 

(M 


=  -  y'  +  fey/v  +  1/v  -  l/a  =  0  , 

(5) 

X  <  X  <  X,  , 

o  —  —  1 

the  prime  indicating  the  differentiation  with  respect  to  x.  The  initial 
conditions  are 

xQ  =  0  ,  v(0)  =  0  ,  y(0)  =  0;  (6) 

the  terminal  conditions  are  not  specified. 


The  quantities  m  and  m  in  (l)  are  bounded,  by  the  inequalities 

m>m.  ,  0<-m<oo,  (7) 

—  man  - 

if  infinite  thrust  is  admitted  as  a  mathematical, convenience.  Such  a  thrust, 
operating  for  an  infinitesimal  time,  produces  arfinite  positive  jump  A  v,  while 
y  and  (<d  -  v)  remain  continuous  in  virtue  of  (5)  and  (4).  In  terms  of  the  new 
variables,  (7)  can  be  written 


\|r,  =  cd 
1  max 


I  I 

y  +  v 


y  -  v  -  x/a  >  0  , 
+  l/a  >  0. 


(8) 
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The  two  unknown  functions  y(x),  v(x)  are  connected  by  a  differential 
constraint  ^  =  0;  the  system  therefore  has  one  degree  of  freedom,  which  can 
be  realized  physically  by  a  choice  of  an  arbitrary  v(x),  ideally  regulated 
by  a  servo-mechanism  controlling  the  flow  rate  m.  Functions  y(x),  v(x)  will 
be  admissible  if  y,  v,  y*  satisfy  the  contraints  (5)  and  (8)  with  the  initial 
conditions  (6),  and  if  they  are  continuous  except  at  corners,  where  y*  and  v 
may  be  discontinuous  with  Av  >  0.  In  the  class  of  admissible  functions  we 

seek  v(x)  that  minimizes  -x  . 
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The  problem  is  thus  identified  with  the  Problem  of  Mayer  in  the 
Calculus  of  Variations,  complicated  by  the  presence  of  algebraic  and 
differential  inequality  constraints. 


3.  THE  AUXILIARY  PROBLEM 

The  differential  constraint  i|r2  >  assuring  the  monotonicity  of  the  mass, 
admits  subarcs  on  which  i|r2  =  0  while  ^  >  0;  i.e.,  the  "burning"  regime  may 
be  interrupted  by  the  insertion  of  "coasting"  subarcs.  In  order  to  avoid  such 
complications  let  us  consider  an  auxiliary  problem  characterized  by  the  absence 
of  the  constraint  \jfg  >  0.  While  such  a  formulation,  used  by  Tslen  et  al., 
automatically  eliminates  the  aforesaid  complication,  it  creates  another  one 
by  admitting  <  0  and  Av  <  0.  Of  course,  negative  fuel  consumption  is  a 
physical  absurdity!  The  resulting  solution  would  not  be  of  physical  interest, 
were  it  not  for  the  curious  fact  that  such  an  occurrence  is  precluded  in  certain 
practical  cases.  Indeed,  if  the  constraints  are  satisfied  anyway  in  the  form 
\|fg  >  0,  Av  >  0,  it  is  clear  that  the  auxiliary  and  the  actual  problems  have  the 
same  solutions.  In  particular,  that  such  is  the  situation  in  both  cases 
treated  here  will  be  shown  in  Theorems  1  and  2  of  section  13.  In  terms  of  the 
quantities  a  and  v,  the  two  cases  can  be  respectively  characterized  by: 

1)  a  is  sufficiently  small; 

2 

2)  C^v  is  convex. 

Accordingly,  we  shall  attack  the  Auxiliary  Problem,  which  is  in  the 
standard  form  of  the  Problem  of  Optimum  Control: 
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"We  seek  a  function  u(x)  satisfying 


i  *  -  y  +  e(x,y,u)  =  o  , 


(9) 


(10) 


o  -  -  1 

subject  to  the  boundary  conditions  and  the  inequalities: 
xq  =  a  ,  y(xQ)  =  b, 
i|(x1,  y^))  =  0  , 
t(x,y,u)  >  0  , 

and  minimizing  some  prescribed  function 

G(x1,y(x1)).  "  (11) 

Here  y,  u,  £  ,  \|r  are  vectors  of  n,  m,  p,  r  components  respectively,  with 
p  <  n  +  1.  In  our  problem  n  =  m  =  r  =  1,  p=0;G=-  x^,  u  =  v,  and 

g  =  (fey  +  l)/v  -  l/a  , 

*  =  “max  "y  “V  ~x/a  =  *1  • 

Since  v  has  disappeared  from  the  problem,  v  has  assumed  the  role  of  a  "control" 
variable,  which  enters  g(x,y,v)  non-linearly.  That  the  problem  is  non-singular 
is  shown  in  section  8;  the  solution  is  obtained  in  sections  5  -  11  by  the 
application  of  the  Necessary  Conditions  I  -  IV  and  the  Fundamental  Sufficiency 
Condition  of  Weierstrass.  The  first  one  is  the  Multiplier  Rule,  comprising 
the  Euler,  the  Transversality,  and  the  Comer  Conditions,  treated  respectively 
in  sections  5,  6,  and  J. 

The  existence  and  the  character  of  the  solution  intimately  depend  on  the 
nature  of  the  drag  coefficient  C^(v) ,  which  is  the  subject  of  the  next  section. 


(12) 


4.  SOME  PROPERTIES  OF  THE  DRAG 

We  shall  assume  the  usual  positiveness  and  the  continuity  of  Cp(v) ,  the 


monotonicity  of  the  drag, 

_d 

dv 


4  (v2)  > 0  » 


(15) 
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and  the  asymptotic  expansions 

2 

Cn  =  A  +  A,v  +  A „v  +  ...  as  v  -»0  , 

D  o  1  2 

CD  =  Bq  +  Bj/v  +  Bg/v2  +  ...  as  v  -»oo  , 

^  0  ^  B^  ^  0  y  i  =  0,  1 f  . .  oo . 

Then  the  logarithmic  derivatives  k  and  k*  defined  by 
k  =  d  log  Cp/d  log  v  ,  k*  =  dk/d  log  v 
have  the  properties: 

k(0)  =  k  (oo)  =  k'  (0)  =  k‘ (oo)  =  0  , 
k(0+)  >  0  ,  k(oo-)  <  0  , 
k'(0+)  >0,  k'  (oo-)  >0  , 
k  +  2  >  0  . 


(l*) 


(15) 


(16) 


FIG.  I  LOGARITHMIC  DERIVATIVES 
k  (v)  AND  k'(v) 
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Furthermore,  let  Cp(v)  have  a  single  maximum  at,  say  aQ.  Then  k  has  a 
maximum  at  a^,  a  zero  at  aQ,  and  a  minimum  at  a^,  while  k*  has  zeros  at 
a1  and  Eg.  It  follows,  in  view  of  (16),  that 

(ao-v)k  >  0  ,  (v  -  a^v  -  a2)k*  >  0  ,  (17) 

0  <  <  aQ  <  a2  <  oo  . 

In  the  analysis,  the  function  f(v),  defined  in  (4),  and  the  derived 
functions  H(v),  h(v),  defined  herewith,  will  be  extremely  useful: 

f  =  CD(v)v2eV/p  >  0  , 

H  =  vfy-  f  ,  (18) 

h  =  H  -  afy  =  (v  -  a)  fy  -  f  , 

with  literal  subscripts  denoting  the  argument  of  differentiation. 

In  terms  of  k  and  k' ,  these  functions  and  their  derivatives  can  be 
exhibited  as  follows: 

fy  =  (f/v)  (2  +  v  +  k)  >  0, 

f  -  (f/v2)  [(2  +  v  +  k)  (1  +  v  +  k)  +  v  +  k']  , 

H  =  f (1  +  v  +  k)  ,  (19) 

\  =  Vfw  » 

h  =  f  [(l  -  a/v)  (2  +  v  +  k)  -  l]  , 

h  =  (v  -  a)f 
v  '  w 

Special  properties  of  these  functions,  obtained  with  the  aid  of  (l6),  are 
tabulated  on  the  following  page: 
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f (0)  =  fv  (o)  =  H(0)  =  Hv(0)  =  h(0)  =  0  , 
f(oo)  =  fy(oo)  =  f^oo)  =  H(oo)  =  h(ao)  =  oo  , 

\  ^(0)  =  CD(0)/p  >  0  ,  (20) 

hy(0)  =  -  a  iw(0)  <  0  , 

p 

2  f  =  f  (0)v  +  ....  as  v  -*  0  . 

w'  '  ’ 

2 

2  H  =  f  (0)v  +  ....  as  v  -> 0  . 

5.  THE  EULER  EQUATIONS 

Since  the  Iagrangian  function  of  the  Auxiliary  Problem  is  F  =  X^  +  p.i|r, 
the  extremals  must  satisfy  the  equations 

y'=  g(x,y,u), 

X'+  Xgy  +  nty  =  0,  (21) 

Xgu  +  U»u  =  0, 
nt  =  0  ,  I|f  >  0  , 

where  X(x),  p.(x)  are  Lagrange  multipliers.  The  subtitution  from  (12)  into  the 
Euler  equations  (21.2)  and  (21.5)  now  yields,  in  view  of  (l8.2), 

X*  +  Xfe  y/v  -  n  =  0,  (22) 

(x/v2)  (Hey  -  1)  -  p  =  0, 

leading  to 

X  =  x(0)  exp  f  £(H  -  vf)  ey  -  lH  dx/v2  ,  (25) 

•'O 

The  use  of  the  "switching  function"  n(x)  permits  simultaneous 
consideration  of  subarcs  lying  in  the  region  >  0,  where  p  =  0,  and  of 
subarcs  lying  in  the  boundary  y1  =  0,  where  p  ^  0.  Three  regimes  are 
distinguished,  designated  by  I,  B,  and  C: 

I.  Impulsive  thrust,  Av  f  0, 

B.  "Burning",  +  >  0,  n  =  0, 

C.  "Coasting",  y^  =  0,  n  /  0. 
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An  extremal  is  compounded  of  a  B- subarc,  with  impulsive  thrusts  I  occurring 
at  a  finite  number  of  points,  and  a  C- subarc  appearing  in  the  terminal  stage 
only. 

During  the  burning  stage  n=0,  and  the  "optimality"  condition  (22.2)  yields 

eyH(v)  =  1  .  (24) 

That  a  solution  v(y)  of  (24)  exists  follows  from  (20),  which  gives  the  range 
of  H  as  (0,  oo  )j  that  this  solution  is  unique  will  be  shown  in  section  9»  with 
the  aid  of  Condition  II.  Several  conclusions  can  now  be  drawn.  First,  (8.1) 
implies  y  <  ooj  then  from  (24)  and  (l8)  there  follows 

H>  0,  v  i  0,  (25) 

and  therefore  v  >  0.  Since  the  initial  value  v  (0)  =  0  violates  the  require¬ 
ment  (25),  the  burning  stage  must  be  preceded  by  an  impulsive  launching  with 
a  velocity  v  that  satisfies  (24)  with  the  initial  condition  y(0)  =  0;  i.e., 

H(vo)  =  1  .  (26) 

The  initial  discontinuity  is  thus  specified  by  v  (0)  =>  0,  v  (0)  =  v  ,  and 
A  a)  =  A  v.  Of  historical  interest  is  the  value  of  the  gravity-drag  ratio  mg/D, 
which  equals  l/fey  in  our  symbols,  and  is  optimized  by  (24)  into  H/f,  or 

mg/D  =  1  +  v  +  k  .  (27) 

The  solution  of  the  Euler  equations  is  obtained  by  the  differentiation 
of  (24)  with  respect  to  x,  yielding 

-  Hy'  =  HyV  ,  (28) 

followed  by  the  substitution  from  (24)  into  (5),  which  in  view  of  (l8)  now 
becomes 

-  a  Hy'  =  h  .  (29) 

Then,  from  (28)  and  (29), 

V  =  h/a  Hy  ,  (50) 

and,  provided  h  /  0,  v(x)  is  obtained  by  the  inversion  of  the  quadrature 

x/a  =  \  dH/h  (51) 

=  X(v)  -  X(vQ)  , 
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where 


(52) 


X(v)  = 


j; 


dH/h 


defines  a  "rocket  function"  dependent  only  on  the  form  of  CD(v)  and  on  the 
value  of  the  parameter  a.  The  equation  of  the  extremal  subarc  now  appears 
in  the  parametric  form  x  =  x(v),  y  =  y(v),  in  consequence  of  (24)  and  (31). 
The  special  case  h  =  0  is  solved  in  section  12. 

During  the  coasting  stage  ijr  =  0,  p  ^  0,  and  (5)  becomes 

0 


w*  +  f(v)  exp  (-v  -x/a  +  to  )  +  1 

IucLX 


(55) 


with  the  initial  conditions  corresponding  to  the  "burnout",  i.e.  the  solution 
of  the  equation  ^  =  0  with  y(v)  and  x(v)  furnished  by  (2k)  and  (31). 

It  is  to  be  noted  that  on  the  C-subarc,  \|r^  =  0,  the  number  of  degrees 
of  freedom  being  zero,  no  variations  are  admitted.  The  implications  are  that 
the  usual  requirement  p  <  0  does  not  hold,  and  that  the  necessary  and 
sufficient  conditions  are  trivially  satisfied. 


6.  THE  TRANSVERSAL3TY  CONDITION 
In  the  control  problem  of  section  3  the  relation 


|(Gx  +  X  .  g)  dx  +  (Gy  -X)  •  dyj  ^  =  0 


(54) 


must  hold  for  all  dx  and  dy  satisfying  the  differentiated  equation  $  =  0,  the 
dot  placed  between  vectors  indicating  their  inner  product.  In  the  Auxiliary 
Problem  n  =  1,  G  =  -  x,  ,  and  p  =  0,  so  that  dx  and  dy  are  arbitrary,  and  (34) 
reduces  to 


1  +  Xg  =  0  , 
X  =  0 


(55) 


at  x  =  x^.  Three  conclusions  can  be  drawn.  First,  since  both  X  and  p  cannot 
vaniBh  simultaneously,  p(x^)  ^  0,  so  that  t^(x^)  =  0.  Second,  noting  that 
|s(x1)l  =  00  from  (35) »  and  recalling  that  y<oo,a>0,  f<oo  ,we  deduce 
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from  (12)  that  v(x1)  =  0.  Of  course,  both  conclusions  are  physically  obvious:  x1 
must  be  reached  with  zero  velocity  after  coasting  with  fuel  consumed.  The 
third  conclusion, 

\(xx  -  0)  >  0  (36) 

follows  from  the  observation  that:  1)  As  v  -»0  the  asymptotic  value  of  g  is 
g <vl/v  >  0,  in  view  of  (12),  (20),  and  (25),  and  that  2)  lim  (Xg)  =  1  as  x  ->x1, 
in  view  of  (35).  Now,  since  x(x)  cannot  change  its  sign  in  virtue  of  (25),  the 
inequality  (56)  implies 

X(x)  >  0  for  xQ  <  x  <  x1  ,  (37) 

which  requirement  can  be  satisfied  by  choosing 

1(0)  =  1  .  (38) 

For  future  use ,  we  note  the  following  asymptotic  values  as  v  ->  0 : 

g  ~  l/v  ,  X  ~  v  ,  X*  -  l/v  ,  \i  ~  -  l/v  ,  (59) 

which  can  now  be  obtained  from  (35)>  (22),  and  (20). 

The  existence  and  the  continuity  of  the  multipliers  X(x),  p(x),  required 
by  the  Multiplier  Rule,  is  now  assured  between  corners  of  a  minimizing  arc. 

7.  THE  CORNER  CONDITION 

At  a  "free"  corner,  the  relations 

A  (X.g)  =0  ,  AX  =  0  (40) 

must  hold,  with  A  denoting  a  jump.  Noting  that  in  our  problem,  with  n  =  1, 

(40)  implies  Ag  =  Ay'  =  0,  and  recalling  that  Ay  =  0  in  virtue  of  (5),  we 
deduce  from  (24),  (29),  and  (18)  that 

AH  =  Ah  =  Afy  =  0  (41) 

on  a  B- subarc.  The  definition  of  H  now  implies 

Af/Av  =  fy  ,  Af  =  0  , 

from  which  the  transition  values  v_  and  v+  can  be  determined. 
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The  transition  values  at  the  corner  are  then  determined  as  follows: 

t 

v_  =  v1  ,  V+  =  Vg  if  v1  >  0 


(46) 


v_  =  Vg  ,  v+  =  v1  if  Vg  <  0  • 

It  will  be  shown  in  section  9  that  such  a  Jump  in  velocity  is  required  by  the 
Weierstrass  Condition  II  whenever  the  Corner  Condition  is  satisfied. 

The  results  of  the  last  paragraphs  can  be  easily  generalized  for  any  N, 
with  (44)  replaced  by 


(47) 


there  being  a  velocity  Jump  for  each  one  of  the  N  double  tangents. 

At  the  Junction  of  the  B  and  C- subarcs  the  Corner  Condition  is  satisfied 

with 

Av  =  Ay  =  Ay=  AX,  =  Ap.  =  0  ,  Av  <  0  .  (48) 


8.  THE  HILBERT  CONDITION 


With  n  =  m  =  r  =  1,  the  four  unknown  functions  y(x),  X.(x),  u(x),  *i(x)  are 
related  by  the  four  equations  (21).  The  highest  derivatives  being  (y*  ,  X* ,  u,  \i) , 
the  non-vanishing  of  the  Jacobian  determinant  is  the  Hilbert  Condition 


F 

uu 


t  0  > 


(49) 


or 

|fuu|  /  o  if  t  >  o  .  (50) 

Here  F  is  the  Lagrangian  function,  and  F^u  is  generally  an  m  x  m  matrix.  The 

condition  assures  the  existence  of  the  highest  derivatives  listed  above,  as 

k-2  k 

well  as  their  piecewise  continuity  of  class  C  if  g  and  +  are  of  class  C  , 

and  Is  a  direct  consequence  of  the  Legendre  Condition  III' . 

In  our  problem  (50)  becomes 


kf^/vH  +  0  , 


(51) 
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and,  since  X,  v,  and  H  are  positive  by  (25)  and  (57), 

fw  I*  0  •  (52) 

Provided  this  requirement  is  met  on  the  B- subarc,  Hilbert  Condition  is  satisfied, 
and  since  g  is  analytic  in  our  problem,  the  functions  (y,X,v,|i)  are  analytic 
between  corners.  It  is  noteworthy  that  the  use  of  the  velocity  v  as  a  control 
variable,  in  place  of  the  thrust  vuf  ,  removes  the  apparent  singularity  of  the 
original  problem. 


9.  CONDITIONS  OF  LEGENDRE  AND  WEIERSTRASS 

The  necessary  conditions  III  and  II,  modified  by  the  inclusion  of  the 
control  variable  u  among  the  set  of  slope  functions,  can  be  written  for  the 
one  dimensional  case,  n  =  m  =  1,  as 

[X6uu  y»  u)  +  ^uu1  5u2  ^  °» 

-  (53) 

E  s  X  [g(x,  y,  u)  -  g(x,  y,  u)  ]  >  0, 

for  all  (x,  y,  u,  X,  n)  belonging  to  the  minimizing  extremal,  and  for  all  u  /  u 
and  satisfying  $  =  0.  In  our  problem  (55),  with  the  aid  of  (12)  and  (24),  reduces 
to 


Xfw/vH  >  0  , 


(x/vH)  [f  -  f  -  (v-v)fv)  >  0, 

where  f  =  f(v),  and  finally,  since  v,  H,  and  X  are  positive,  to  the 
requirements  that 


(54) 


f 

w 


>  0 


> 


f  -  f  >  (v  -  v)fv 


hold  on  every  B- subarc. 


(55) 
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In  the  language  of  geometry,  (55)  implies  that  v  must  be  restricted  to 
the  domain  where  the  tangent  to  the  curve  f(v)  lies  entirely  below  the  curve. 
For  drag  of  type  1,  (55)  is  automatically  satisfied;  for  drag  of  type  2,  where 
(44)  and  (45)  hold,  (55)  is  equivalent  to  the  requirement 

(v  -  vx)  (v  -  v2)  >  0  ,  (56) 


where  v^  and  Vg  are  the  points  of  contact  of  the  double  tangent.  The  exclusion 
of  the  interval  (v1>  Vg)  from  the  B-subarc  then  demands  that  a  corner  occur 
when  v  reaches  the  values  v^  or  Vg,  as  described  in  (46).  Conversely,  the 
occurrence  of  such  a  corner  satisfies  the  requirement  E  >  0,  the  equality  holding 
only  at  corners  for  v  =  v^,  v  =  Vg,  and  conversely.  Consequently,  E  >  0  except 
on  a  set  of  measure  of  zero,  so  that  the  Weierstrass  Condition  holds  in  its 
strengthened  form  II ' . 


A  fortiori,  the  strengthened  Legendre-Clebsch  Condition  nf  ;  i.e., 


f  >  0 

W 


(57) 


also  holds,  from  which  two  consequences  follow.  First,  (57)  establishes  the 
Hilbert  Condition  (52);  second,  with  the  aid  of  (19.4)  it  implies  that  ^  >  0. 

We  conclude  that  H(v)  is  monotonic  in  the  domain  defined  by  (56),  and  has  an 
inverse  H~  ,  thus  assuring  the  uniqueness  of  the  solution  v(y)  of  the  equation  (24). 
In  view  of  this  fact,  it  is  convenient  to  replace  H  and  h  in  all  the  equations 
referring  to  the  B-subarc  by  H*  and  h*  (see  Fig.  5)  defined  by 

H*  =  H,  h  =  h*  if  (v  -  V;L) 

H*  =  H(v^)  ,  h*  =  h  (v1)  if  (v  -  vx) 

In  the  future,  if  no  confusion  results,  the  asterisks  of 


(v  -  v  )  >  0  , 


(58) 

H  and  h  will  be  dropped. 


(v  -  v2)  <  0  . 


10.  THE  JACOBI  CONDITION 


The  Second  variation  in  the  control  problem  defined  by  (9)  -  (ll)  with 
n=m=r=l,  p  =  0  can  be  written 


d2J  = 


f 

J  x 


F  )  dx2  +  2F  dydx  +  d2G 

y  J  ^ 


(59) 


(Fyy  6y  +  2Fyv6y5v  +  F^Sv2)  dx  , 
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where  F  is  defined  by  F  =  +  (h|r  ,  i  =  -  y  +  g(x,  y,  v)  ,  (it  =  0,  t  >  0.  The 

necessary  Jacobi  Condition  is  that 

d2J  >  0  (60) 


must  hold  on  a  minimizing  extremal  for  all  dx,  dy,  and  for  all  By,  6v 
satisfying  the  differentiated  equations  ^  »  0,  (it  *  0. 


Observe  that  F  =  -  X*  from  (21.2),  and  that  our  problem  possesses  the 
y  2 

following  special  features:  l)  G  =  -  x^,  so  that  d  G  =  0;  2)  Fx  =  -  \i/a  in 
view  of  (12);  3)  u  =  v;  and  F  =  since  ^  is  linear  in  v;  4)  By  and  Bv 

must  satisfy 


-  By'  +  gyBy  +  g^v  =  0  , 
(i  (By  +  Bv)  =  0 


(61) 


Furthermore,  on  the  B-subarc,  xq  <  x  <  £,  where  £  is  the  "burnout"  point,  observe 
that  (i=0,  gy  =  0  by  (21.3),  and  the  initial  condition  y(o)  =  0  implies  By(0)  =  0, 
and  hence  the  solution  of  (6l)  is  By  =  0,  Bv  arbitrary.  On  the  other  hand,  on 
the  C-subarc,  §  <  x  <  x^,  (t  ^  0,  and  the  continuity  of  By  in  (6l)  implies 
6y(|  +  0)  =  8y(g  -  0)  =  0.  Hence  the  solution  of  (6l)  is  By  =  Bv  =  0,  in 
agreement  with  the  remark  in  the  last  paragraph  of  section  5.  Now,  since 
By  =  0  everywhere,  dy  =  y'dx,  and  (59)  becomes 

d2  J  =  Jj-|i/a  -  X*g)  dx2[j  +  f  Xg^&v2  •  (62) 

J  x 

♦  1  o 

Finally,  at  x  =  x1  (39)  yields  (i<v-  l/v  and  X*  g~-  l/v  as  v  ->  0;  on  the 
B-subarc  Xg_  >  0  by  (53. l)  and  (57).  Since  a  and  v  are  positive,  we  conclude 

vV  |  g 

that  the  Jacobi  Condition  is  satisfied  in  its  strengthened  form  IV  ,  d  J  >  0. 

An  immediate  consequence  is  that  in  our  problem,  with  n  =  1,  IV1  assures 
the  existence  of  a  field.  Let  a  family  of  extremals  y(x,0)  be  defined  by 
(9)  -  (ll)  with  the  initial  condition  (10.1)  replaced  by 

xq  =  a  ,  y(xQ)  =  b  +  8  ,  (63) 

where  8  is  a  family  parameter.  That  the  region  bounded  by  x  *  a,  $  (x,y)  =  0 
is,  indeed,  a  field  follows  from  the  following  considerations:  1)  The 
extension  of  IV*  to  0  /  0  is  trivial;  2)  if  assures  the  simple  covering  of 
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the  region;  i.e.,  the  existence  of  the  function  0  (x,y) ,  and  hence  of  the  slope 
functions  u(x,y)  and  multipliers  X(x,y),  p(x,y);  3)  in  a  one-dimensional 
problem,  the  Euler  equations  suffice  to  assure  that  the  Hilbert  integral  is 
independent  of  the  path. 


11.  THE  SUFFICIENCY  CONDITION 

We  resort  to  the  following  variant  of  the  Fundamental  Sufficiency  Condition 
of  Weierstrass,  proved  in  Appendix: 

"Let  a  family  y(x,0)  of  extremals  of  a  control  problem  be  generated  by 
the  initial  conditions  (63),  involving  0  as  a  parameter.  If  this  family 
constitutes  a  field,  and  if  each  extremal  of  the  field  satisfies  I  and  II*  with 
the  appropriate  initial  conditions,  then  the  extremal  for  0=0  yields  an 
absolute  minimum  of  the  control  problem." 

Note  that  Conditions  I  and  if  have  been  established  in  sections  5-9  for 
9=0,  and  that  their  extension  to  the  family  defined  by  (63)  is  trivial. 
Furthermore,  the  existence  of  a  field  has  been  proved  in  section  10.  We 
conclude  that  the  hypothesis  of  the  theorem  is  satisfied,  and  that  our  extremal 
therefore  yields  an  absolute  minimum  of  the  Auxiliary  Problem. 

12.  THE  STEADY  SLATES  OF  MOTION 

Aside  from  their  intrinsic  interest,  the  lenmas  of  this  section  are 
required  in  the  proof  of  the  Basic  Theorems  of  section  13. 

* 

Lemma  1:  "The  function  h  (v)  has  one  and  only  one  positive  zero,  7." 

The  proof  proceeds  from  (18),  (19),  (20),  (57) >  (58) •  Two  cases  are 
distinguished: 

case  1.  a  is  outside  (v^,  Vg)  (See  Fig.  3) 

Then  the  relations  h  =  (v-a)  fy  -  f,  hy  =  (v-a)  f  ,  f  >  0  imply  that  h* 
has  one  and  only  one  stationary  point, 

min  h*  =  h(a)  =  -  f (a)  <  0.  (64) 
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*  # 

Since  the  minimum  is  negative,  the  relations  h  (0)  =0,  h  (oo) 

-X* 

continuity  of  h  imply  the  conclusion  of  the  lemma,  with 
a  <  y  <  v1  or  v2  <  a  <  y  , 

(v  -  y)h  >  0  . 


oo,  and  the 


(65) 


case  2.  a  is  inside  (v^,  Vg) 

Then  h*,  stationary  on  the  interval  (v^,  Vg),  attains  there  a  minimum, 
min  h*  =  h(v^) .  That  the  minimum  is  again  negative  is  implied  by  h(o)  =  0, 
h(a)  <  0,  and  0  <  v1  <  a;  finally,  h(v1)  =  h(vg)  <  0  and  if  (oo)  =  oo  imply 
the  conclusion  of  the  lemma,  with 

v2  <  a  <  v2  <  y  ,  (66) 

(v  -  y)  h*  >  0  . 

Lemma  2;  "On  the  burning  subarc  of  an  extremal,  the  acceleration  of 
the  rocket  cannot  chanf  2  its  sign." 

The  proof  proceeds  from  (50)  and  (19),  leading  to 

v'  =  h/a^  ,  hv  =  (1  -  a/v)  ^  , 

(67) 

x 

h  =  hQ  exp  J  (l/a  -  l/v)  dx  . 

Noting  that  a  >  0,  and  that  >  0  by  (19),  (25),  and  (57),  we  conclude  that 

9  *  * 

sgn  v  =  sgn  h  =  sgn  hQ  .  (68) 


The  Corollary,  "hQ  =  0  implies  h(x)  =  0  and  v  =  7"  follows  immediately. 
Three  types  of  trajectory  are  thus  distinguished: 


a)  vQ  <  7  , 

h(x)  <  0  , 

v*  <  0, 

deceleration  ; 

b)  vQ  >  7  , 

h(x)  >  0  , 

1 

v  >  0, 

acceleration ; 

O 

<J 

O 

II 

>• 

h(x)  =  0  , 

< 

in 

steady  state; 

In  case  c)  the  solution 

(24)  and  (51), 

of  the  Euler  equations,  must  be  replaced 

by  y  «  0,  v  =  7. 
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13.  THE  BASIC  THEOREMS 


Having  constructed  the  solution  of  the  Auxiliary  Problem,  we  shall  show 
that  under  the  assumptions  of  Theorems  1  or  2  it  satisfies  the  constraints 
>  0  and  Av  >  0  on  the  B- subarc. 

Theorem  1:  "2a  <  min  (a^  v1)  implies:  1)  a  <  y  <  2a,  2)  > 

3)  Av  >  0,  with  the  last  two  relations  holding  on  the  B- subarc  of  an  extremal 
of  the  Auxiliary  Problem. 


To  prove  1),  observe  that  from  (18.3) 


h(a)  =  -  f(a)  , 

h(2a)  =  f (2a)  [  a  +  \  k  (2a)] , 


(69) 


and  that  the  hypothesis  and  (17)  imply 

0  <  2a  <  a,  <  a  , 

1  o 

(70) 

k(2a)  >  0  . 

Then  (69)  and  f  >  0  imply  h(a)  <  0  and  h(2a)  <  0.  Furthermore,  from  the 
hypothesis,  a  <  2a  <  vx  <  vg,  so  that  h  =  h*  on  (a,  2a),  in  view  of  (58).  The 
conclusion  is  now  implied  by  the  continuity  of  h*  and  by  1, 

To  prove  2),  observe  that  y2,  defined  by  (8.1)  as 

t»». 

V2  =  =  y  +  v  +  l/a  ,  (71) 

can  be  exhibited  as  a  function  of  v,  with  the  aid  of  (29),  (30),  (l8),  (19), 
in  two  alternate  forms: 


♦2(v)  =  h/aHy  +  fy/H 


(72) 


=  H/o^  +  (fv/vHy)[2  +  k  +  (v  +  k')/(l  +  v  +  k)]  , 


The  positiveness  of  v,  f,  H,  fy,  Hy,  2  +  k,  1  +  v  +  k  is  assured  by  (16),  (18), 
(19),  (20),  and  (57).  There  are  two  possibilities:  Either  h  >  0  or  h  <  0,  If 
h  >  0,  then  ^  >  0  in  the  first  line  of  (72),  On  the  other  hand,  if  h  <  0,  then 
v  <  7  in  (65);  the  previous  conclusion,  7  <  2a,  and  the  hypothesis,  2a  <  a1, 
imply  v  <  a,j  then  (17)  implies  k  >  0,  and  finally  >  0  in  the  last  line 
of  (72). 


2k 


To  prove  3),  recall  that  v1  >  2a  >  7,  and  that  by  Leona  2 ,  v  >  7  implies 
h  >  0,  v*>  0,  and  conversely.  Then  note,  with  the  aid  of  (4l),  that 
h(v1)  =  h(v2)  >  0.  Therefore  v*  >  0  if  v  =  v1  or  v  =  Vg,  thus  excluding  the 
possibility  Av  <  0  in  the  second  line  of  (46). 

For  rough  practical  purposes,  the  hypothesis  of  the  theorem  may  be  replaced 

by 

KM  >  4  ,  (73) 

where  M  is  the  jet  Mach  number,  and  K  is  the  ratio  of  the  atmospheric  specific 
heats.  To  derive  this  result  observe  that:  l)  a  =  gl/c2  =1/  K  M2  by  the  law 
of  perfect  gas  and  the  formula  for  the  sonic  velocity,  2)  the  constants  a^  and  v^, 
which  are  the  values  of  v  at  the  maximum  of  k(v)  and  at  the  first  point  of 
contact  of  the  double  tangent  to  f(v),  respectively,  lie  in  the  sonic  region, 

3)  the  sonic  velocity  corresponds  to  aQ~  l/M,  and  both  a.^  and  v^  are  generally 
sufficiently  near  aQ  to  justify  the  inequalities  a1A0  >  l/2  ,  v1/aQ  >  l/2. 

Thus,  (73)  implies  2a  <  a^and  2a  <  v^ 

For  the  Earth,  with  c  ~  2000  m/s,  g~9.8m/s  ,  i~8000  m/s,  K  =  1.4,  we 
calculate 


a  ~  0.02,  M  6, 
KM  a/ 8.4, 


(74) 


concluding  that  the  hypothesis  of  the  theorem  is  satisfied  for  terrestrial 
rockets. 

The  vacuum  case,  p  =  0,  solved  by  Mlele,  corresponds  to  a  =  0  and  is, 
therefore,  a  subcase  of  the  theorem.  On  the  other  hand,  the  constant-density 
atmosphere,  p  =  const.,  corresponds  to  a  =  oo  and  hence  lies  outside  the  scope 
of  Theorem  1.  Indeed,  Leitman  succeeded  in  solving  this  case  only  by 
invoking  the  quadratic  law  of  drag,  CD  =  const.,  bringing  the  problem  within 
the  scope  of  Theorem  2. 

2 

Theorem  2:  "If  G^v  is  convex,  then  ^g  >  0  and  Av  >  0  on  the  B-subarc  of 
an  extremal  of  the  Auxiliary  Problem." 

In  the  proof,  note  that  the  hypothesis,  in  view  of  (15),  implies 

(k  +  1)  (k  +  2)  +  k’  >  0,  (75) 
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Mid,  consequently,  t2  >  0  111  the  last  line  of  (72).  Furthermore,  observe  that 
(75),  (16. 4),  and  (19.2)  imply  the  convexity  of  f(v).  Hence  the  drag  belongs 
to  type  1  of  section  7,  with  no  comers  on  the  B- subarc,  and  with  Av  =  0. 

The  cases  of  quadratic  law  of  drag  and,  the  more  general,  power  law  of 
drag,  also  treated  by  Miele,  appear  as  subcases  of  the  theorem. 

14.  SUMMARY 

Under  the  assumptions  of  Theorems  1  or  2  the  solution  is  characterized 
by  the  structural  formula 

(IB)N  +  1C;  (76) 

i.e,,  the  burning  stage  B,  preceded  by  an  impulsive  launching  I,  contains  N 
additional  impulsive  thrusts,  N  being  the  number  of  double  tangents  of  the 
curve  f(v),  and  is  followed  by  the  coasting  stage  without  fuel.  The  solution 
therefore  Includes  as  a  special  case  the  results  of  Tsien  and  Evans,  where 
N  =  0.  An  absolute  minimum  has  been  established  with  the  aid  of  the  second 
variation  and  a  variant  of  the  Sufficiency  Principle  that  is  particularly  useful 
in  problems  of  optimum  control. 

Theorem  1  applies  to  terrestrial  launching  and  any  drag  function  with 
some  very  general  properties  listed  in  section  4;  Theorem  2  covers  extra¬ 
terrestrial  launching  but  is  restricted  to  a  fairly  common  class  of  drags  that 
includes  all  the  cases  previously  treated  in  the  literature. 

BORIS  GARFINKEL 
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APPENDIX 


To  prove  the  Sufficiency  Condition  stated  in  section  11,  define  v(x,y)  by 
w  =  G(x,y)  +  I*  , 


* 

I 


J{u- 


g(x,y,u)  +  nt]  dx  -  X.dy 


}. 


(77) 


* 

where  I  is  the  Hilbert  integral  of  the  control  problem  and  u,  A,  n  belong  to 
the  field.  With  A  denoting  an  increment,  observe  that:  1)  Aw  =  0  on  a 
closed  path;  2)  Aw  =  0  on  a  boundary  subarc  in  $  =  0,  in  virtue  of  the 
Trans versality  Condition  (34)  and  ^\|r  =  0;  3)  Aw  =  AG  on  an  extremal  of  the 
field,  in  view  of  y  =  g(x,y,u),  where  u(x,y)  is  the  "slope  function". 


Next  let  0  and  1  denote  the  end-points  of  the  extremal  for  ©  =  0,  and 
let  CQ2  be  any  admissible  arc  connecting  0  and  a  terminal  point  2  lying  in 
<£=  0.  Then  there  follows  from  the  properties  of  w  listed  above  that 

wQ1  =  G(l)  -  0(0)  , 


(78) 


wq2  =  0(2)  -  0(0)  +  I  (CQ2)  , 


+  w. 


leading  to 

0(1)  -  0(2)  =  I*  (co2)  .  (79) 

Finally,  note  that,  in  view  of  (53)  and  |4r  =  0,  the  expression  for  I  in  (76) 
can  be  also  exhibited  as 


I*  =  j*  X.  (g  -  g)  dx  =  -  J*  E  dx  , 

S  —  g(x,  y,  u)  ;  u  ^  u, 

f  * 

and  that  II  implies  E  >  0  on  CQ2,  and  hence  I  <  0  in  (80)  and  (79). 
conclusion 

0(1)  <  0(2) 


(80) 


The 

(81) 


follows  immediately. 
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